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We present a new method to calculate the total ion-ion interaction potential in terms of building 
blocks which we refer to as "single-particle interaction potentials". This allows also to compose 
the separate contributions from neutrons and protons to the interaction potentials. The method is 
applied to nuclear collisions via the use of time-dependent Hartree-Fock theory. 
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The knowledge of interaction potentials between com- 
posite many-body systems is of fundamental importance 
for predicting the outcome of reactions involving such 
entities. These are often modeled using empirical func- 
tions that depend on macroscopic variables, such as the 
distance between the two centers, treating the composite 
objects as structureless particles and ignoring the mi- 
croscopic origins of these potentials (e.g., Lennard- Jones 
potential for rare-gas atoms nucleon-nucleon poten- 
tials 0). On the other hand, microscopic approaches 
try to obtain such potentials by including the interac- 
tions of the constituents as the building-blocks and cal- 
culating the whole potential as a function of some set of 
macroscopic variables (e.g., Born-Oppenheimer approxi- 
mation for molecules deformed shell models for nu- 
clei However, the total potential obtained this way 
does not reveal the contribution coming from individual 
single-particle states. Quoting Schrodinger, best knowl- 
edge of a whole does not necessarily include best possible 
knowledge of the parts. The entanglement of these con- 
tributions may yield further insight as to the interaction 
of many-body systems via a representative potential. 

Most common microscopic approaches for calculating 
interaction potentials usually employ the adiabatic or 
sudden approximations for the relative motion of the in- 
teracting systems @. Recently, we have introduced a 
new microscopic approach for the calculation of ion-ion 
potentials for nuclear collisions. This method is based 
on the time-dependent Hartree-Fock (TDHF) descrip- 
tion of nuclear collisions coupled 6] with a constraint 
on the local density at the given instant of time. In- 
teraction potentials obtained using the so-called density 
constrained TDHF (DC-TDHF) method have been 
successful in describing fusion cross-sections for a num- 
ber of systems |8|. The new approach we will describe 
below is generally suitable for microscopic calculations 
provided Koopmans' theorem Q is applicable. 

The microscopic approaches based on the mean-field 
approximation often use the energy difference between 
the combined system and the asymptotic energies of the 
individual systems to calculate the ion-ion interaction po- 
tential as 

V{R) = E Al+M {R)-Ef i -Ef i , (1) 



where Ea 1 +a 2 {R) 1S the total energy obtained for the 
combined system as a function of the distance R be- 
tween the two ions, and E A °^ are the individual bind- 
ing energies of the two systems calculated using the 
same interaction. The binding energies are often cal- 
culated either by integrating the energy density over 
all space E = J d 3 rT-l(r) or via Koopmans' expression 
E = ^ J2 a w a{ta+£a) , where w a denotes the occupation 
probability of the single-particle state a. This mixture of 
half kinetic energies t a and half single-particle energies 
e Q applies to Hartree-Fock calculations employing purely 
two-body forces. For mean-field calculations based on 
density-dependent effective forces |10| or energy- density 
functionals fill ] a rearrangement term needs to be added 
to the above expression [12| ■ This can be disentangled to 
single-particle energies and rearrangement term as 




(2) 

where v is the density-dependent part of the interaction. 
We can deduce a generalized single-particle sum from 
that by decomposing the second p in the rearrangement 
term into single-particle densities as p — Y] a p a . This 
yields 

E = ]^^w a (t a +e a + e a ) , (3) 

a 

where e a is the single-particle rearrangement energy aris- 
ing from the above modification. All of the binding en- 
ergies comprising the ion-ion interaction potential via 
Eq. ([l} can be expressed using the above form. In terms 
of these expressions the ion-ion interaction potential can 
be represented as a sum over single-particle interaction 
potentials as 

V(R) = J2 w ^a(R) , (4) 

a 

with 

v a (R) = ~ [(t Q +e a + ~e a )(R) - (t° + e° + e°)] , (5) 

where the quantities with a superscript (0) denote the 
asymptotic values. This decomposition is sensible if the 
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single particle energies are good representatives of the 
corresponding separation energies according to Koop- 
mans' theorem Q and if the rearrangement corrections 
remain small. Such a situation is given in the application 
example discussed below, namely nuclear mean field cal- 
culations using Skyrme forces. In that case, the proton 
and neutron interaction potentials are obtained simply 
by summing the single-particle potentials v a (R) over all 
protons, or neutrons respectively. For electronic energy- 
density functionals one has to employ a self-interaction 
correction to restore Koopmans' theorem approximately 

[HQ 




FIG. 1. (color online) Interaction potential for a head-on 
collision of 16 0+ 16 at E c .m. = 12 MeV. The black curve 
is the total ion-ion interaction potential. The red and blue 
curves show the contribution to the total potential coming 
from protons and neutrons, respectively. Also shown is the 
point Coulomb potential. 
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FIG. 2. (color online) Contribution of neutron single-particles 
to the interaction potential for a head-on collision of 16 0+ 16 
at Ec.m. = 12 MeV. The sum of all these potentials give the 
neutron curve in Fig. [Tj The arrow marked Rb indicates the 
barrier top of the total potential shown in Fig. [T] 



In practice, we have used the Skyrme SLy4 effective 
nuclear interaction for our calculations. The rear- 
rangement terms described above can be written as; 



For the application of the above concepts we have used 
the DC-TDHF approach Q. In this approach TDHF 
time-evolution takes place with no restrictions. At cer- 
tain times during the evolution the instantaneous densi- 
ties p p / n (r,t) are used to perform a static Hartree-Fock 
minimization while holding the neutron and proton den- 
sities constrained to be this instantaneous TDHF densi- 
ties [Hj]. In essence, this introduces the concept of an 
adiabatic reference state for a given TDHF state. The 
difference between these two energies represents the in- 
ternal energy. The adiabatic reference state is the one 
obtained via the density constraint calculation, which is 
the Slater determinant with lowest energy for the given 
density with vanishing current and approximates the col- 
lective potential energy after the subtraction of the static 
binding energies as shown in Eq. ([!} . All of the dynami- 
cal features included in TDHF are naturally included in 
the DC-TDHF calculations. 
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e 2 / d 3 r p a pp , a <E {protons} 



d 3 r p a X n , P (r) 



where we have defined p a = ^* (r)^> a (r) , and A(r) is the 
coordinate dependent Lagrange multiplier for the density 
constraint [15| . The subscripts n,p in the rearrangement 
term for the density constraint indicates the choice cor- 
responding to the isospin content of the index a. The 
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FIG. 3. (color online) Same as Figure [2l but for protons. 

summation over a exactly reproduces the total rearrange- 
ment energies obtained by direct integration. 

The first collision studied is 16 0+ 16 at E c . m . = 
12 MeV. In Fig. Q] we show the total interaction potential 
(black curve) as well as the contributions coming from 
neutrons and protons to this total potential. The to- 
tal interaction potential is numerically identical to the 
one obtained in standard DC-TDHF calculations using 
the integral of the energy density to obtain the terms in 
Eq. ([1]). As a reference we also show the point Coulomb 
interaction. As can be seen from the figure the outer 
part of the potential barrier is primarily determined by 
the interaction potential between the protons of the two 
nuclei while the neutron potential is essentially zero. In 
other words the neutron-proton interaction does not in- 
fluence the outer barrier region. While inside the barrier 
neutrons provide all of the attraction and the proton po- 
tential remains positive for all R values. 

In Fig.[5]we show the individual neutron single-particle 
potentials for the 16 0+ 16 system. In spite of using 
three-dimensional Cartesian coordinates, we are able to 
calculate some quantum numbers when the solution pos- 
sesses a good symmetry. In this case spherical initial 
16 nuclei allow the use of spherical quantum labels at 
the asymptotic position of the two nuclei. On the other 
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FIG. 4. (color online) Interaction potential for a head-on 
collision of the 16 0+ 24 at E c . m . = 12 MeV. 



hand, when the nuclei overlap the system becomes ap- 
proximately axially symmetric thus allowing the compu- 
tation of parity and the z-component of total angular 
momentum. The reassignment of parities and the be- 
havior of the states as they evolve towards smaller R 
values is noteworthy. For smallest R values the filling of 
the levels resembles the one for a single-centered 32 par- 
ticle system. The trend seen in the actual single-particle 
energies is relatively easy to understand; the two oxygen 
nuclei are initially far apart and are brought together. 
While they are far apart and isolated the nuclei have 
identical energy levels. However, as the separation be- 
tween the two nuclei becomes smaller the single-particle 
states begin to overlap. The Pauli principle dictates that 
no two nucleons in an interacting system may have the 
same quantum state. Therefore, each (doubled) energy 
level of the isolated nuclei splits into two orbitals, one 
lower in energy than the original level and one higher. 
This can be most easily seen for the two initial lsi/2 
states. However, in addition to energy splitting the par- 
ity of one of these states also changes from positive to 
negative since no more than two neutrons can be in the 
lsi/2 state (each state is originally occupied by two neu- 
trons, w a = 2). In the language of molecular physics the 
states that attain a lower energy than their asymptotic 
value are referred to as the bonding states, whereas the 
states that evolve to a higher energy are the anti-bonding 
states Q. We observe that not all neutron states are 
bonding states but to the contrary about half are actu- 
ally anti-bonding in character for small R values. One 
of the P3/2 states (green curve) make a transition from 
anti-bonding to bonding for smaller R values. Another 
manifestation of the bonding and anti-bonding states is 
their spatial localization. Bonding states are localized in 
the interior part of the combined system as opposed to 
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anti-bonding states that are more spread out, as can been 
seen from the single-particle moments. Figure [3] shows 
the contribution of proton single-particle states to the to- 
tal ion-ion potential. The trend of the proton states are 
essentially the same as the neutron states with the excep- 
tion of the rise from the zero potential line (dotted lines) 
due to the presence of the Coulomb interaction. The 
sum of all these proton single-particle potentials repro- 
duces the point Coulomb potential for the two incoming 
nuclei as can been observed in Fig. [TJ We have also per- 
formed calculations for head-on collisions of 16 0+ 24 at 
E c . m . = 12 MeV, and for 40 Ca+ 40 Ca at E c . m . = 55 MeV. 
In Fig. [4] we show the neutron and proton potentials for 
the 16 0+ 24 system. We note that in comparison to the 




12 14 
R(fm) 



FIG. 5. (color online) Same as Figure [4] but for a head-on 
collision of the 40 Ca+ 40 Ca at E c . m . = 55 MeV. 
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O+ O system the potential minimum is considerably 



lower due to the presence of eight extra neutrons. The 
surplus bending from the extra neutrons also affect the 
proton potential, bringing it down to negative values for 
small 7?. We also note that the neutron potential starts 
its dip to negative values at larger R values than the 
16 0+ 16 case thus bringing the total barrier maximum 
to a lower value. Figure [5] shows the same quantities for 
the 40 Ca+ 40 Ca system. The behavior of the proton po- 
tential is significantly different in this case. Inside the 
barrier the proton potential is essentially constant and it 
only rises as the minimum value of R is reached. This 
is probably due to the stronger Coulomb potential which 
counter-weights the nuclear binding. While the neutron 
and total potentials appear intuitively as expected, the 
proton potential behaves in an unexpected manner. This 
is further evidence that the total potential may be mask- 
ing some interesting features of its building blocks. 

We have introduced a general approach for the cal- 
culation of single-particle interaction potentials as the 
building blocks of the total interaction potential for two 
fragments of finite fermion systems. The formulation is 
general provided the Koopmans' theorem is applicable 
to the underlying energy functional used for the many- 
body calculations. The formalism is applied to calculate 
ion-ion potentials for nuclear reactions using TDHF for 
the time-evolution of the nuclear collision together with 
the density constraint formalism to find the correspond- 
ing adiabatic reference state. We show the contribution 
of single-particle potentials to the total potential for the 
16 0+ 16 case. We identify repulsive and attractive con- 
tributions as bonding and anti-bonding states. Perhaps 
the more obvious aspect is the identification of neutron 
and proton contributions to the total potential where we 
see clearly the subtle interplay between Coulomb repul- 
sion and nuclear attraction for the protons while neutrons 
are always the dominant contributors to binding. 
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